Exact result for nonreciprocity in one-dimensional wave transmission 
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For sound waves impinging on a one-dimensional medium, we show that nonlinearity can lead to 
nonreciprocal transmission, without dissipation or broken time reversal invariance. Placing quasi- 
monochromatic filters at the ends of the nonlinear medium, nonreciprocity can be obtained without 
the generation of higher harmonics outside the medium. Remarkably, in this configuration the 
nonreciprocity is found to be proportional to the net energy flow when monochromatic sources of 
equal strength (at the fllter frequency) are simultaneously turned on at both ends. This result is 
conjectured to be general for one dimensional scattering. It is also shown that although simultaneous 
monochromatic sources lead to net energy flow, with sources of small but non-zero bandwidth there 
is no net energy transport, in accordance with the second law of thermodynamics. 

PACS numbers: 43.25.-|-y, 46.40.Cd 
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The reciprocity theorem has a long history in acoustics 
and optics. For the case of a linear medium with time 
reversal invariance, the theorem can be proved HQ: for 
one dimensional systems, it amounts to the transmission 
coefficient being the same when waves are incident from 
the left or from the right. In the absence of time re- 
versal invariance, reciprocity is no longer necessary. In 
optics, this can be achieved by intrinsic magnetization 
in the scatterer or with an external magnetic field 
In acoustics, nonreciprocity can be caused by — and be 
used to detect — the motion of objects, as in acoustic 
tomography Q. 

For nonlinear media, however, even if time reversal 
invariance is not broken, nonreciprocity is possible 0. 
Photonic structures with diode like behavior have been 
proposed, where the effect of the nonlinearity is strength- 
ened by the existence of a bandgap 0, Q . Such passive 
diode like behavior can be useful in the field of optical 
communications . 

In this paper, we examine nonreciprocity in one- 
dimensional nonlinear media for longitudinal waves such 
as sound. As with light, we find that nonlinearity is suf- 
ficient to result in nonreciprocity, even for dissipationless 
systems without broken time reversal invariance. We go 
on to consider a different configuration: when monochro- 
matic filters are placed at the two ends of the nonlin- 
ear medium. This confines the higher harmonics to the 
medium, so that the reflected and transmitted wave are 
at the same frequency as the incident wave. Thus there 
is no 'contamination' from higher harmonics outside the 
nonlinear medium. This setup, apart from possible ad- 
vantages from a communications perspective, allows us 
to examine constraints from the second law of thermo- 
dynamics. 

With filters, we obtain the unexpected result that the 
nonreciprocity is now proportional to the net energy 
transport from one side of the system to the other when 
two monochromatic sources of equal strength are simul- 
taneously connected to the two ends of the device. We 



conjecture that this proportionality is general for any 
scattering process with two input and two output chan- 
nels (at the same frequency) that is invariant under time- 
translation and time-reversal and is perturbatively acces- 
sible (explained later in this paper). 

We consider a system that can be modelled as two adja- 
cent layers, in each of which longitudinal waves propagate 
in accordance with a nonlinear wave equation. Outside 
the system, both to the left and the right, the linear wave 
equation is satisfied. Thus we have 



(1) 



where y is the displacement of the wave, and ni,Bi,^i 
vary from region to region. Inside the scatterer, the two 
layers have parameters (ni, Bi, fii) and (rt2, ^2, ^2)- Out- 
side the scatterer, n = B = 1 and fi ~ 0. The scat- 
terer is taken to cover the region — 1 < a; < 1, with the 
boundary between the two layers at x = 0. The form of 
Eq.(^) retains the leading order nonlinearity in the elas- 
ticity of the medium; the energy density of the wave is 
^n'^y^ + ^{dxy)^+It{dxy)^/3. At the three boundaries be- 
tween the four regions, y and Bd^y + ii{dxy)^ (the force 
exerted on the boundary from the two regions it sepa- 
rates) are continuous. We also consider an alternative to 
Eq.O 
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B.dly + ti,dx{dxyf 
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which is slightly easier to work with, but which has an 
accidental y —y symmetry. Eqs.flJ and (0) are in the 
class of Fermi Pasta Ulam (FPU) wave equations |3] ■ 

For both Eqs.QJ and ((2Jl, we first use perturbation the- 
ory to obtain an analytical solution. The incoming wave 
amplitudes from the left and right are ai and 02 respec- 
tively. Eqs.l^l and ||2Jl can then be solved to linear order 
in ai,2, and then iteratively to successive higher orders 
in perturbation theory. For the case without filters, one 
imposes the requirement that all frequency components 
of the solution are purely outgoing outside the scatter- 
ing medium, except for the component at frequency uj 
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whose incoming part is specified. For the case with fil- 
ters, except for the component at frequency uj which is 
unaffected by the filters, all other components are con- 
fined to the scattering medium and have zero amplitude 
at a; — ±1. These conditions are sufficient to solve Eqs.flJ 
and 10), order by order. 

The equations were solved to third order using 
Mathematica"'"^ , with = 1 and various specific val- 
ues chosen for ni_27^i.2 and w. This third order solu- 
tion yields the leading 0(|a|^) correction to the outgoing 
power to the left (or to the right) for Eq.(Q. On the other 
hand, a similar third order solution to Eq.|(21) yields the 
outgoing power to 0(|a|^), with two nonlinear contribu- 
tions to the component at frequency lo. 

For the case with filters, the outgoing wave is entirely 
at frequency lo. Expressing the outgoing amplitudes 61^2 
as an expansion in powers of ai_2 and a* 2, time trans- 
lational invariance requires that each term in the expan- 
sion should have one extra power of the unconjugated 
variables as compared to the conjugated ones. Thus 

i>i{ai,a2) = M^jUj + Nijkia*akai + The outgoing 

power to the left is It is possible to verify for both 

Eq.Q and Eq.© that |&i(0, a)p - |62(a, 0)p is not equal 
to zero, demonstrating nonreciprocity in the transmission 
coefficient. (Since the system is nondissipative, this is 
equivalent to nonreciprocity in the reflection coefficient.) 
With /^i.2 — 1 and various values of ni^2, ^1,2, and u all 
- 0(1), the coefficient of \a\'^ in |6i(0, a)|2 -'|62(a, 0)|2 is 
0(1). In units where /i = 1, the amplitude of the incom- 
ing wave will be small, so this is a weak effect. However, 
it should be possible to enhance the effect by construct- 
ing more complicated structures, for instance acoustic 
analogs of Ref. 6]. For the case without filters, one has 
to consider the outgoing power at different harmonics 
separately, but nonreciprocity is still found. 

With filters, one might consider a monochromatic 
source as a blackbody (white noise) source from which 
only waves of one frequency are allowed to escape If 
one connects a blackbody at each end of the nonlinear 
scatterer, with both blackbodies at the same tempera- 
ture, there should be no net flow of energy from one side 
to another. With the fllters, this would seem to be equiv- 
alent to choosing ai^2 to be equal in magnitude, but with 
a random relative phase. Surprisingly, it is possible to 
verify through the perturbation expansion of the previ- 
ous paragraph that if 01^2 are indeed of equal strength, 
the phase averaged outgoing power is not the same on 
both sides of the scatterer. In fact, for all choices of the 
parameters we have tried, we have verified that 

|6i(0,aV2)|2- |52(aV2,0)|2 
= 2(|6i(a, ae'^)p)^ - {Ma, ae*^)P)^. (3) 

This has been verified to third order for both Eqs.QJ and 
©, i.e. the 0{\a\'^) terms for Eq.©, and the 0{\a\^) 
and 0(|a|^) terms for Eq.Q. We have been unable to 
find any result resembling Eq.Q for the case without 
filters, either including or excluding higher harmonics in 



the outgoing wave. 

Since our analytical treatment is only perturbative, we 
turn to numerical simulations. The nonlinear medium 
is modelled by a chain of N particles with anharmonic 
springs connecting them. Thus if are the displacements 
of the particles from their equilibrium positions, 

"^^^/^ = -dy^[V{yi - + V{yi+i - y^)] (4) 

for all the particles except the first and the last one, with 

V{y)^ly' + ly'. (5) 

The first and the last particle must be coupled to the ex- 
ternal environment. This coupling is through incoming 
and outgoing waves, with — as in any scattering prob- 
lem — the incoming waves specified and the outgoing 
waves determined by the scattering medium. Beyond 
the left boundary of the medium, the external waves can 
be expressed as fi{x — vt) + fo{x + vt). The force ex- 
erted by these waves on the boundary of the medium is 
proportional to —dx[fi + Jo]- By continuity, the veloc- 
ity of the boundary is equal to the velocity just outside 
the scattering medium, which is dt[fi + /o]- From the 
form of fi and /o, it is easy to see that —vdx[fi + /J = 
—dt [fi + /o] + 29t/i. Thus the external force acting on the 
boundary is a sum of a term proportional to the velocity 
of the boundary, and a term specified by the incoming 
waves. For monochromatic waves, we have 

rrnyi = -miujlyi - V'{yi - j/2) - + M cos{ujt) 
rriNyN = -friNUilyN - V'{yN ~ Vn-i) 

—KyN + A2 cos{ujt + (j)) (6) 

where <j) is the relative phase between the incoming waves 
from the left and the right. Thus the coupling to the 
external environment is seen as an effective damping and 
forcing term in the equation of motion for the first and 
last particle. 

The flrst term on the right hand side of Eqs.@ makes 
these particles act as filters ii lo = ljq and muv are very 
large. Due to the nonlinearity of the medium the in- 
coming waves at frequency wq produce a response at all 
multiples of loq. If the excitations are resolved into fre- 
quency components, for the component at wq the left 
hand side of Eqs.© cancels the first term on the right 
hand side. The forcing and effective damping term from 
the external environment must balance the dyV{y) term 
from the interior, as they would have if the terminal parti- 
cles had been missing. On the other hand, at any higher 
harmonic, mi^N^uj'^ — loq) diverges in the mi^pf — s- 00 
limit, so that yi^Ninujo) for 71 ^ ±1. Thus for the 
component at loq the terminal particles are transparent , 
whereas for higher harmonics the terminal particles act 
as fixed boundaries for toi^at — > 00, confining the higher 
harmonics to the nonlinear medium. 

Eqs.Q and © together with Eq.© were numerically 
simulated for a chain of 4 -I- 2 particles, in units where 
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FIG. 1: Numerical results for a 4 + 2 particle chain with 
asymmetric mass-distribution. The particle masses are 100, 
1.7, 1.4, 1.9, 1.3 and 100. The left-to-right current Jio = 
162(1,0)1^ is plotted as a function of the nonlinearity e. The 
nonreciprocity, Jio — Joi, is shown for 7 different values of e, 
and compared with the phase averaged current Jn with both 
sources on (right vertical scale). The two quantities are equal 
in the perturbative state, but not in the non-perturbative 
state. 



K — ujQ — I and Ai_2 — 1- Various values of r7i2 • • • 
were used; the results shown in Figure 1 are representa- 
tive. The nonlinearity parameter e was varied, which is 
equivalent to a fixed e and varying Ai^2- The chain was 
started at e = 0, allowed to reach steady state, and then 
e was increased slowly till approximately e — 0.4, after 
which it was decreased slowly to zero. The system under- 
goes a transition as e is increased, from a periodic state 
with frequency lo — loq ^ I to a noisy state. The transi- 
tion between the two is first order, with an accompanying 
hysteresis loop, but for sufficiently small or large e only 
one state is seen. The existence of two states is similar to 
that in Ref. for optics. Perturbation theory, which is 
'connected' to e = and only allows for harmonics of coq, 
cannot access the noisy state. As seen in Figure 1, the 
perturbatively accessible periodic state satisfies Eq.Q, 
whereas the nonperturbative state does not [T^ . 

In the noisy state, preliminary results when the inci- 
dent wave is entirely from the left (on mi) show broad 
peaks in the transmitted power at w w OAujq and ~ 
O.lSwo. As e is increased further, there is another transi- 
tion between e « 0.35 and 0.45, with a jump in the av- 
erage transmitted power and a broadband component to 
the power spectrum [l^ . The jumps from the perturba- 
tive state to the first noisy state and thence to the second 
noisy state are at different values of e when the incident 
wave comes from the left instead of the right. A detailed 



dynamical analysis would be required to characterize the 
various noisy states and the transitions between them. 
However, this is not the focus of this paper. 

The non-zero right hand side of Eq. ^ might seem to 
contradict the second law of thermodynamics. If two 
blackbody sound sources at the same temperature were 
connected at the ends, no net energy flow would be pos- 
sible. The filters would only allow waves at frequency 
luq to enter or exit the system, seemingly equivalent to 
monochromatic sources. However, as can be seen from 
our numerical implementation, any filter has a non-zero 
(albeit arbitrarily small) bandwidth. For a nonlinear 
medium, the different frequency channels interact with 
each other. Thus even with filters, a blackbody and 
monochromatic source are not strictly equivalent 'l4'|. 
This may seem like a quibble, but from the discussion 
before Eqs.® it is clear that blackbody sources at the 
ends would correspond to white noise being applied to 
the terminal particles (which is then filtered by them). 
Eqs.Q) and © are then generalized Langevin equations 
(with damping and noise only in Eq.©), which can be 
rigorously proved to reach thermal equilibrium |l5j . In 
view of our explicit results for monochromatic radiation, 
and the Langevin description for blackbody sources, we 
must conclude that narrow and zero bandwidth filters 
are not equivalent beyond linear order |l6j| . The situa- 
tion here is different from the one considered in Ref. ^3 j 
where non-equilibrium energy sources were used; since 
the sources had to be maintained out of equilibrium, sec- 
ond law arguments were inapplicable there. 

We note in passing that Eq.Q) is the standard FPU 
system .9], which is difficult to equilibrate ^3|> but the 
open boundaries in Eq.® seem to be sufficient to equi- 
librate the system with thermal (blackbody) sources. 

We return to the possible basis of Eq.©- With fil- 
ters, the scatterer can be viewed as generating a map- 
ping from the two complex input amplitudes to the two 
complex output amplitudes. This mapping has to satisfy 
the properties that i) since the system is nondissipative, 
+ = |aiP + |a2p ii) from time translation invari- 
ance, if ai^2 — > 01,26*" then bi^2 — > ^1.26*" hi) from time 
reversal invariance, if 01^2 ~^ ^1 2 then 61^2 ~^ 2 i^) the 
mapping is perturbatively accessible from the zero am- 
plitude limit. We conjecture that these constraints may 
be sufficient to yield Eq.Q. This would imply that the 
equation is valid for any one dimensional two-channel 
scattering problem that satisfies the conditions above. 
In view of our numerical results, the fourth condition 
is essential; mappings that violate Eq.® can in fact be 
constructed without it p9|. 
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